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■ Abstract. Cosmic-rays are ubiquitous, but their origins are surpris- 

ingly difficult to understand. A review is presented of some of the basic 
Q ' issues common to cosmic particle accelerators and arguments leading to 

Q , the likely importance of diffusive shock acceleration as a general explana- 

tion. The basic theory of diffusive shock acceleration is outlined, followed 
by a discussion of some of the key issues that still prevent us from a full 
understanding of its outcomes. Some recent insights are mentioned at 
the end that may help direct ultimate resolution of our uncertainties. 
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! The inherent difficulty in understanding the acceleration of cosmic-rays (CRs) 

may not immediately be obvious. At the most basic level we must presumably 
identify an electric field capable of producing particles of very high energy. That 
sounds straightforward in fast moving plasmas. For galactic and especially for 
ultra-high energy CRs, the energies involved are so large that the possibilities 
are very limited. When we consider, in addition, the energy distribution of 
53 . the CRs, as well as their composition, rate of production and other details, 

however, the task of modeling their production and propagation becomes very 
sophisticated. In this talk I will deal mostly with a few of the more common 
^ I and basic issues as they apply to baryonic galactic CRs below the "knee" , which 

■ we can conveniently take to be ~ 10^'^GeV/nucleus. Several speakers at this 

meeting have admirably addressed many of the special issues relevant to other 
aspects of the broader problem. 

There is now broad consensus that galactic CRs are accelerated mostly from 
the interstellar medium (ISM) at supernova remnant blast waves by the diffu- 
sive shock acceleration (DSA) process. Beyond that simple statement, however, 
significant differences of opinion quickly surface on almost every detail. Despite 
decades of concerted and highly productive effort, this is not yet a solved prob- 
lem, either physically nor astrophysically. I will now briefly outline some of the 
arguments pointing us to the consensus viewpoint for the basic scenario, then 
follow with a brief outline of a few of the issues that continue to hinder our 
efforts to solve the problem fully. While these complicating issues seem to be 
major barriers to a comprehensive understanding, there are hints that when all 
the pieces of DSA theory are in place together, a robust and possibly simple 
product may result. For additional DSA insights I direct readers to the accom- 
panying discussion by Kang (2001), which focuses on empirical evidence for DSA 
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as well as some of important numerical and technical issues and how they are 
being addressed. 



2. Background Issues 

Ultimately, CR acceleration comes through an electric field; however, the most 
convenient descriptions may not show this explicitly for a given process. The 
electric fields are most likely inductive, through large scale motions, although 
they may be directly applied through stimulated plasma waves. In any case 
we can express their effective magnitude as f ~ PaB, where Pa is the relevant 
speed in the accelerator and B is the strength of the local magnetic field. For 
an accelerator of length scale, Ra: 

we can use this to constrain the necessary 

magnetic field as 

where is the required particle energy in units of lO^GeV and Z is the charge 
on the particle. This simple constraint can also be derived through a number of 
different conceptual approaches with modest variations in the numerical constant 
and some variation in the interpretation of Pa- Examples include using equation 
[5] for the time needed for DSA to produce the required energies, constraining the 
diffusive length scale of a particle to be smaller than the size of the accelerator, 
or even just requiring the particle gyroradius to be smaller than the size of the 
accelerator. 

Figure 1 illustrates the result of equation 1 for Ej = 1 and two values of 
Pa in a form made popular by Hillas (1984). In the figure I have indicated 
rough conventional model properties for a small sample of astrophysical objects. 
The only galactic objects known that may be able to satisfy the constraint 
are supernova remnant shocks ("SNRs"), winds from O and B stars, pulsar 
magnetospheres and possibly compact accreting binaries. 

Pulsars can be excluded as the principal source of galactic CRs by consider- 
ing CR composition. As described in detail by other speakers at this workshop, 
the CRs below the "knee" roughly mirror the composition of the sun and the 
ISM (e.g., Seo 2001; Wiebcl-Sooth, Bicrmann & Meyer 1998). There arc impor- 
tant differences, including at the isotopic level, that provide vital clues about 
the details of the source plasma and CR propagation history (e.g., Meyer, Drury 
Sz Ellison 1997). But, to lowest order this information tells us that the source 
material is almost surely the ISM, with perhaps some small admixture of locally 
processed stellar material. 

Beginning from that point most models of galactic CR acceleration have 
focused on SNRs based on the total energy input required. That can be esti- 
mated by considering the rate at which CR accelerators must replace CRs that 
diffuse from the galaxy. Isotopic ratios fix the characteristic escape time near 
100 MeV to be ~ lO^yrs (Council 1998). Taking the observed local flux at the 
solar system and assuming this is also the average for the galaxy this leads to an 
average required CR power input near 10^-*^ erg/sec (e.g., Drury, Markiewicz & 
Volk 1989; Fields et al. 2000). If we take for comparison the galactic supernova 
rate to be ^yr~^, and kinetic energy yield per event to be lO^^erg, the avail- 
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Figure 1. Plot after Hillas (1984) showing a general constraint be- 
tween the size of a particle accelerator, its characteristic magnetic field 
and the maximum energy possible in the accelerator. The diagonal 
lines are lower bounds identified in relation [1] for a proton energy 
of lO^^eV and two different characteristic velocities in the accelerator. 
Rough properties are estimated for several types of accelerators. 
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able power in SNRs is rough 10^^ erg/sec, which is sufficient, but not by a huge 
factor. Thus, SNR-based models must be at least moderately efficient to supply 
the needed power. No other known galactic source comes closer than an order 
of magnitude to this power supply, explaining why most models have involved 



The CR flux energy spectrum below the knee measured at earth is a power- 
law after correction for solar modulation. 



where (j) is the energy flux, and f{p) is the phase space distribution function 
(e.g., Seo 2001; Wiebel-Sooth et al. 1998). These scalings are strictly valid only 
when the CRs are relativistic. In addition, the flux, is very nearly isotropic, 
consequent to the diffusive propagation of CRs through the ISM. Propagation 
models also generally lead to a steepening of the spectrum with respect to its 
form at the source, by an increment ~ 0.5 — 0.6 in the slope (e.g., DuVernois, 
Simpson & Thayer 1996), reflecting an energy dependence to the apparent CR 
escape rate. Thus, CR source models usually aim to explain a power-law distri- 
bution function, f{p) a p^'', with g « 4.1 — 4.2. The fact that the simple steady 
state test particle DBA theory predicts a power-law f{p) with q 4+ when 
shocks are strong is one of the primary reasons that model for CR acceleration 
has attracted so much attention over the past two decades. I will, in fact, limit 
my remaining discussion to issues associated with this process. 

3. An Outline of Diffusive Shock Acceleration Theory 

There are a variety of approaches to understanding the physics underlying DSA, 
since the microphysics is complex and depends on what one assumes about 
such details as the structure and orientation of the local magnetic field. All 
approaches depend on a small fraction of nonthermal particles becoming trapped 
by scattering around a shock front, so that they may tap into the energy flow 
through the shock for extended times, but with a finite probability of escaping 
in a given time interval. It is remarkable that all these approaches give virtually 
the same answer to the first approximation, at least so long as feedback on the 
plasma flow can be neglected and various pathologies are avoided. For in-depth 
discussions of the theory there are a number of flne reviews (e.g., Berezhko & 
Krymskii 1988; Blandford & Eichler 1986; Drury 1983; Jones & Ellison 1991; 
Malkov & Drury 2000). Here I will outline only some basic properties of the 
theory as it applies to quasi-parallel shocks, since that is one of the simplest 
situations. 

The DSA theory depends on an almost isotropic particle distribution and 
particle speeds large compared to the bulk flow speed, u. Spatial gradients lim- 
ited by the assumption of diffusive propagation with respect to local scattering 
centers lead to a transport equation of the diffusion-convection type (e.g., Parker 



SNRs. 



(I){E) oc p^f{p) oc p"^-^ oc E' 



1-2.7 
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1965; Skilling 1975a), 



^+u-Vf = -^p^ (V • n) + V • (nVf) + Q, 



(3) 
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where the first term on the right accounts for adiabatic compression, the second 
spatial diffusion and Q is a generic source term that can represent injection 
or escape, for example. The diffusion coefficient is given by k = |Aw, where 
A is the scattering length, and v is the particle speed. The flow velocity is u. 
In this context the scattering is commonly assumed to involve resonant Alfven 
waves, whence one can estimate from quasi-linear theory (e.g., Skilling 1975b) 
a scattering length A = (^Vg, with Vg the particle gyroradius and 

Here Pb is the total magnetic pressure and kPyjk is the pressure (energy density) 
in waves satisfying resonance, which can be "sharpened" to be expressed as 
kp = LOcm, or kvg = 1, where lOc is the nonrelativistic cyclotron frequency for the 
particle species under discussion. The limiting value, C, = 1 {\ = Vg) leads to so- 
called Bohm diffusion. For an oblique magnetic field at a plane shock this same 
formalism applies with the substitutions k ^ Ky, then k = Kycos^^ -|- K.^siv?0, 
with 9 the angle between the magnetic field and the shock normal and kj_ = 
K,\\/{l+C,'^) (e.g., Jokipii 1987). Then k refers to diffusion along the shock normal, 
while K|| and kj^ describe diffusion along and perpendicular to the local mean 
magnetic field. 

For a parallel shock [0 = 0) we can imagine the acceleration as a first 

order Fermi process, with particles successively being scattered across the shock 
from opposite sides of a converging flow. The mean fractional momentum gain 

between successive downstream returns is ^ ^^^"^ , where Au = U2 — ui is 
the velocity change across the shock. The probability of downstream escape 
by advection following each downstream return is simply P^sc ~ 'iu2/v. Since 
the average distance a particle diffuses on either side of the shock before being 
returned is Xd^^ = ki, 2/^^1,2, the mean time between crossings is tsc = x^/v. 
Using this one can estimate the mean time for a particle to be accelerated from 
pi to p2 as (e.g., Lagage & Cesarsky 1983) 
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— ^ (factor) X = (factor) x td{p2), (5) 
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where the arrow represents a trend to the Bohm limit, and the "factor" is 20 for 
a limiting strong gas shock, if ki = K2, and p2 » pi- Requiring ta < R/u\, we 
recover relation [1] multiplied by a factor with = ^ in the Bohm limit. 
Note that Xd/ui = t^, called the diffusion time for the CRs, is the time scale 
over which an isotropic population of CRs will be advected across the diffusive 
"precursor" formed ahead of a shock. This is also the average time for CRs to 
diffuse a length Xd = n/ui. 

For a plane gas shock the steady solution to equation [3] is a power-law 

with 

_ 3ui _ 3r 4M2 
^~ U1-U2 ~7^^ M2+3' ^ ' 

where r = ^ = ^ is the compression ratio of the shock, M is the shock Mach 
number, and the arrow corresponds to a 7 = | gas. Then g — > 4 in the strong 
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shock limit. A simple computation shows that q = 3+PescX so kinematically 
the spectrum reflects the match between the rate of particle acceleration and 
escape. This solution neglects any backreaction of CRs on the shock structure, 
so constitutes a "test particle" solution to the problem. 

That simple, limiting solution, apparently independent of any microphysical 
details, and naturally leading very close to the expected source slope for galactic 
CRs was one of the key insights that raised the community's consciousness about 
DSA in the late 1970s (Axford, Leer & Skadron 1977; Krymskii 1977; Bell 1978; 
Blandford &: Ostriker 1978). The other was a realization that the Alfven waves 
needed to isotropize the CRs would be generated by the CRs themselves as they 
streamed into the oncoming upstream plasma. The postshock plasma is expected 
to be turbulent (e.g.. Quest 1988), including waves advected from upstream, so 
scattering in that region seemed assured, as well. Quasilinear theory provides 
an estimate of the growth time for resonant Alfven waves ahead of the shock, 
which depends on V/. Using the fact that the CRs will diffuse upstream a 
characteristic length, x^, so that V/ oc f{p)/xa, we can estimate this time to be 

xd Pu, , Ma Pb 

VAPcijp) C Pcijp) 

Here Pc{p) is the pressure in resonant CRs, va is the upstream Alfven speed, 
and Ma is the Alfvenic Mach number of the shock. Thus, on the surface, this 
time scales with the acceleration time, although the additional factors, which 
are not all known or even constant, complicate the comparison considerably. In 
practice calculations have generally assumed tyj is very short, so that P^ reaches 
an asymptotic limit (commonly given as Bohm diffusion), or that wave dissi- 
pation and growth are in a local equilibrium that produces another preferred 
diffusion coefficient (e.g., Jones 1993). In fact, except for some hybrid plasma 
simulations involving a limited range of particle energies (e.g.. Quest 1988; El- 
lison, et al. 1993) this is not a solved problem. 



4. Some Important Details 

The real beauty of DSA was its apparent simplicity and the robust character of 
the solution. However, virtually as it was introduced, DSA theory exposed some 
potentially very messy details that make the simple behaviors outlined above not 
obviously valid in many applications. One I will mention but not elaborate here 
is the influence of an oblique magnetic field at the shock. Then the magnetic field 
component aligned with the shock face jumps at the shock, and that means there 
is necessarily an electric field also aligned with the shock face, but orthogonal 
to the magnetic field. This leads to "shock drift acceleration" and additional 
possible complications if the particle propagation is not diffusive across field 
lines (e.g., Gieseler et al. 1999). Jokipii (1982) showed, however, that so long 
as the particles diffuse across field lines and the shock is planar, then the basic 
formalism already outlined remains intact even for 9 ~ tt/2. When >> 1, 
however, the length and time scales for highly oblique shocks can be reduced 
from Bohm diffusion values (Jokipii 1987). Another complication can result if 
the magnetic field itself "wanders" or is "braided" (Kirk, Duffy & Gallant 1996), 
since particle motions along the field lines can lead to additional spatial transport 
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that influences the rates of acceleration and escape. While some shocks may well 
be almost perpendicular, with 6 ~ 7r/2, it seems to me if the magnetic field is 
moderately turbulent, that global perpendicular shock behaviors are not likely 
to be prevalent for nonrelativistic shocks. 

Many of the other complications can be summarized by noting that DSA 
is an integral part of collisionless shock formation itself (Eichler 1979). The 
particles we call CRs arc really just a nonthermal tail of the full distribution, 
fip). They are distinguished by their relatively long scattering lengths, and, 
of course, the relatively large individual particle energies, reflecting the absence 
of full thermodynamical equilibrium. The CR particles come from, or "are 
injected" from, the more abundant bulk plasma population. They also exert 
a pressure through their interactions with Alfven waves that can modify the 
flow of the bulk plasma. The strength of those waves depends, in turn, on the 
intensity of the CR streaming, as already pointed out, and on comparative wave 
growth and dissipation rates. Wave dissipation heats the plasma as well. Flow 
modiflcations from a classical gas shock structure resulting from these features 
alter the various terms in the diffusion-convection cqiiation [3] , so a test particle 
solution based on an unmodified shock flow must be reexamined in light of a 
more integrated view of the physics. 

4.1. Dynamical Backreaction 

The possible importance of CR backreaction was quickly recognized from esti- 
mates of the likely CR pressures (Axford et al. 1977; Eichler 1979). That can 
be expressed as 



For f{p) DC this diverges logarithmically as q ^ 4 for strong shocks and as 

P2 ^ oo. Clearly a real shock must at least include a cutoff at finite p2- Such a 
cutoff appears naturally from equation [5] in a shock of finite age, or as a result 
of escape by CRs above some momentum The latter effect may result either as 
a consequence of finite shock extent, or from limitations to the intensities of 
Alfven waves resonant with the highest momentum particles. I will revisit this 
last point later on. 

Using equation [7] we can crudely compare Pc at the shock for a given CR 
number density to the thermal pressure, Pt = UtkhT ~ ritPtVt, where pt = mvt 
represents the (nonrelativistic) momentum of a thermal ion. Taking g = 4 to 
illustrate, and assuming for simplicity that the CRs are all relativistic, we have 



Since all the terms on the right are large except for the fraction of ions in the 

CR populations, ric/nj, it becomes immediately obvious when P2/P1 >> 1 that 
a relatively very small CR population can easily produce a pressure comparable 
to the thermal gas. 

The consequences of backreaction to the shock properties are significant. 
First, the gradient from a finite Pc slows and compresses plasma adiabatically 





Pc Uc pic P2 



(9) 



Pt nt ptvt pi ' 
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before it reaches the classical and discontinuous gas "subshock". This effect 
can be simply estimated by recalling that the CRs are distributed upstream 
in a precursor of characteristic length x^, producing a pressure gradient that 
decelerates the flow as it approaches the shock by an amount 



Aui dPc 1 Xd 1 Pc Xd Pc 



c 



ui dx Pi ui ui Xd Piu\ Piu\ 



(10) 



As Pc at the shock becomes comparable to the dynamical momentum flux into 
the shock, we expect the subshock to become very much weakened, since 

Atti/iti ~ 1, and adiabatic heating in the shock precursor will reduce the Mach 
number of the flow entering the subshock. 

The first and simplest approach to evaluating in detail the dynamics of 
modified CR shocks used the energy moment of the diffusion-convection equation 
(Drury & Volk 1981; Axford, Leer & McKenzie 1982). Defining 

Ec = ATTmc^ p [\f^^ - l) /(P)f^lnp, (H) 

and the closure relation, Pc = (7c — l)-£'c, we have from equation [3] 
BE 

+ u ■ VEc = -7cV • n + V • (< K > VEc) + S, (12) 

provided the limits pi and p2 can be neglected. The term < k > is a mean 
diffusion coefficient weighted by momentum and f{p), and S is an integral form 
of the source term, Q, in equation [3]. In equation [11] I have expressed p in units 
of mc. When merged with Euler's equations for gas dynamics, this approach 
is commonly termed a "two-fluid" dynamical model, since the CRs are treated 
as a massless, diffusive fluid coupled to the bulk flow. Backreaction on the 
bulk plasma is included through the pressure gradient terms in the bulk flow 
momentum and energy equations. The pondermotive force of the Alfven waves, 
as well as their energy dissipation may also be readily included (Achterberg 1982; 
McKenzie &: Volk 1982). The two-fluid approach is somewhat controversial, 
mostly because the closure parameters, 7c and < n > are not known a priori, and 
because of some pathological steady state solutions for strong shocks in the limit 
7c = |. Nonetheless, when properly used, it is an effective and computationally 
efficient method to establish basic features of modified CR shocks (Kang &; Jones 
1995). 

The first two-fluid computations showed that it was possible for most of 
the momentum flux through a shock to be converted into CR pressure, for 
example, with Pc amplified over the precursor length, Xd- As many as three 
steady solutions were identifled for strong shocks from given upstream conditions 
or particle injection rates, when 7^, ss | (Drury Sz Volk 1981). The solutions 
differ substantially in the "efficiency" of conversion by the shock of momentum 
influx, piul, into Pc, and represent a bifurcation phenomenon with respect to 
the supply of seed particles. The pathological steady state solutions involved 
finite postshock Pc from zero upstream Pc, and completely smoothed shocks, in 
which the gas subshock disappears. Those solutions are the result of assuming 
P2 ^ 00 and cannot be reached practically from time dependent solutions, or 
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when p2 is finite. Nonetheless, as I will outline below from other considerations, 
the discovery from two-fluid models is correct that shocks may either be highly 
efficient or not very efficient in accelerating CRs in ways that depend sensitively 
on the supply of seed CRs. Two- fluid calculations also confirmed that the basic 
evolutionary timescale for shock modification is ta = Xd/ui (Drury &; Falle 
1986), and identified the existence of dynamical instabilities derived from the 
long scaled coupling between the bulk plasma and the CRs (Drury Sz Falle 1986; 
Zank, Axford & McKenzie 1990; Ryu, Kang & Jones 1993). 

4.2. Nonlinear Modifications to the CR Spectra 

According to the rightmost expression of equation [6] the CR spectrum should 
steepen when a pressure precursor forms at the shock, so that the Mach number 
of the subshock is reduced by adiabatic heating and deceleration of the inflow. 
That is a somewhat misleading observation, however, since the total compression 
across the structure, including the precursor, is greater than for the gas shock 
alone. Thus, according to the other expressions for q in equation [6] particles 
scattered across the full shock transition, where r > 4 for an initially strong 
shock, would be expected to form into a very hard spectrum with g < 4. The 
relevant interaction length for the CRs is, x^, of course. For Bohm-like diffu- 
sion, with K oc pv, we have Xd{p) oc pv, so particles at relatively low momenta 
respond mostly to the jump across the gas subshock, while the highest momenta 
particles will reach well past the diffusion lengths of those lower momentum par- 
ticles. From relations [9] and [10] it is clear that modest momentum CRs may 
produce significant compression in front of the subshock. That property tends 
to produce concavity in the form of fip), and that, in turn enhances the relative 
importance to Pc of the highest momentum particles over the expression [9] ; that 
is, the divergence of the pressure is faster than logarithmic with p2- Thus, the 
hydrodynamical form of the shock and the form of the CR distribution, f{p), are 
linked in a highly nonlinear manner. This point is crucial to our understand- 
ing of DSA in practice, as pointed out by a number of authors (e.g., Ellison 
&; Eichler 1984; Malkov 1999). The development of these nonlinear features in 
such a modified CR shock is clearly visible in the evolution of the fully nonlinear 
diffusion-convection-equation-based simulation presented in (Kang 2001; Kang 
et al. 2001), for example. 

It is not yet entirely clear what a strongly modified CR shock will look like 
when it is examined in a complete and fully self-consistent way, nor what the 
CR spectrum is, despite a considerable effort put into determining those issues. 
I will return at the end to some recent insights into those questions. Before 
that it is useful to complete our discussion of issues with a few comments on 
two more critical aspects of the problem; namely, the injection of CRs out of 
the bulk plasma, and feedback between the CR acceleration and the evolution 
of the Alfvenic turbulence responsible for moderating the acceleration. 

4.3. The Critical Role of Injection 

Several authors have pointed out that the efficiency of CR acceleration at strong 
shocks depends sensitively on the rate of injection there (e.g., Eichler 1979; 
Berezhko et al. 1995; Malkov 1999). Berezhko et al. (1995) argued for the 
existence of a critical injection rate, above which the shocks are highly efficient. 
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so that Pc at the shock is a large fraction of the momentum flux into the shock. 
Below such a threshold the process becomes much less efficient, so the pattern is 
reminiscent of the original two-fluid results. As part of a study of CR acceleration 
in SNRs Berezhko et al. found a sharp increase in acceleration efficiency in 
high Mach number shocks as the injected proton fraction was increased between 
10~^ and 10~^ of the number flux through the shocks. They arbitrarily fixed 
that number in the models, but the simulations still highlight the issue clearly. 
Malkov (1997a, b) also argued for an injection threshold on the grounds that 
once a shock begins to be modified, so that the compression is increased, the 
highest momentum particles see a larger velocity jump, so are more effectively 
accelerated, thus enhancing Pc- That, in turn, enhances the compression at the 
shock, increasing the acceleration rate, and so on. The process then becomes 
limited by the highest momentum to which CRs can be accelerated. 

It becomes crucial, therefore, to incorporate appropriate injection physics 
in DSA models. As noted, CRs are an extension of the thermal particle pool 
reflecting the absence of full equilibrium. At a shock the majority of ions are 
"thermalized" and unable to re-cross the shock, since the shock thickness is de- 
termined by the characteristic thermalization length of the ions. That process 
is very complex and incomplete in a collisionless plasma (Kennel, Edmiston & 
Hada 1985), however, and some fraction of ions having been only partially ther- 
malized may escape upstream as "seed particles" for DSA. The injection problem 
amounts to determining how that seeding, or "thermal leakage" is controlled. 
Monte Carlo simulations handle the process very simply by assigning a form 
to the scattering law for all the ions that allows a smoothly increasing escape 
probability with increasing momentum; i.e., by setting A oc with Ps > (e.g., 
Baring, Ellison & Jones 1994). While that captures the flavor of the process, it 
does not attempt to include an explicit model for the nonlinear plasma physics 
associated with the thermalization process. Hybrid plasma simulations do that 
in detail, of course, but are not really designed to explore the production of very 
high energy particles that may be accelerated at cosmic shocks. 

Malkov (Malkov & Volk 1995; Malkov 1998) has recently developed a very 
promising analytical model for thermalization and associated injection based 
on the nonlinear trapping of ions in Alfven waves generated by ions escaping 
upstream and amplified through the shock. This model is calibrated against 
hybrid plasma simulations, so contains no free parameters. They find a very 
sharp cutoff in the probability to return upstream. For strong shocks only ions 
with speeds more than roughly 10 times the bulk flow speed away from the shock 
have a finite chance to escape back into the oncoming flow and become seed CRs. 
That is a pretty strong filter, and in nonlinear diffusion-convection simulations 
we carried out recently based on this model (Gieseler, Jones &: Kang 2001), 
the injected proton fraction quickly stabilized around 10~^. For this model that 
result seems fairly robust, with an equilibrium formed between reduced injection 
coming from cooling of postshock gas as seed particles escape, against decreased 
upstream compression and a stronger subshock resulting from reduced Pc if 
injection falls below the equilibrium value. 
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4.4. Alfven Wave Feedback Loops 

By now it should be clear how intricately connected the different elements of the 
nonlinear DSA model are. Once injected the rate at which particles are acceler- 
ated in a parallel shock depends on two things; namely, the velocity profile of the 
bulk flow and the spectrum and intensity of resonant Alfven waves across the 
flow profile. As mentioned at the beginning, amplification of the Alfven waves in 
the precursor is generally attributed to instabilities fed by the CRs themselves 
as they attempt to stream ahead of the shock. While the quasilinear theory of 
wave amplification is wcll-cstablishcd, and was used in deriving equation [7], for 
example, once the wave amplitudes become large, quasi-linear theory is suspect. 
Similarly, wave dissipation is generally attributed to nonlinear Landau damping 
(e.g., Volk, Drury & McKenzie 1984), but again that has not been developed to 
a state that it can be reliably used to give an accurate, fully nonlinear treatment 
of source terms for the Alfvenic turbulence. The wave dissipation is important 
in another way, since it leads to local heating of the bulk plasma, adding to the 
adiabatic heating that already reduces the Mach number at the gas subshock. 
This heating acts as another limiter in the acceleration process. 

Malkov, Diamond & Jones (2001) have also pointed out an important detail 
in the transport of Alfvenic turbulence in modified shocks that strongly infiu- 
ences the maximum particle momentum, p2, that, we will recall, becomes the 
controlling influence in the efficiency of DSA at modified shocks. To see this, 
consider at a given time that the maximum momentum for the CRs is p2 (t) , and 
then follow those particles as they subsequently return from the downstream flow 
moving into the upstream flow with increased momenta, p2 = P2(0 + ^Pi after 
scattering. Being the first to stream into the flow with these momenta, they 
do not encounter significant Alfven wave amplitudes at resonant wavelengths 
{krg 1). Thus, they should easily escape the shock and will not be further 
accelerated. Their streaming will, on the other hand, amplify whatever low level 
Alfven waves are upstream at the resonant wavelength, and those waves will be 
advected towards the shock where they can interact with subsequent CRs trav- 
eling upstream. However, in a strongly modified shock the flow is compressed as 
it approaches the gas subshock, causing the in-flowing Alfven wave to be com- 
pressed, as well. (The Alfven speed will generally decrease.) Therefore, these 
waves will be resonant only with CRs of smaller momentum than p2, and the 
current population of CRs at p will not be scattered until they propagate to 
regions where there has been no compression to the flow. Their rate of return 
is consequently reduced by the preshock compression, so that there may be a 
substantial reduction in P2it) from the one predicted from Bohm diffusion. 

5. Discussion: Resolving these Issues 

The preceding section may leave one with less than full confidence that we will 
soon be able to model fully nonlinear modified CR shocks in a complete and 
self-consistent manner. There are, however, some encouraging developments 
that could lead us into a much clearer understanding. For one, computational 
techniques, as discussed by Kang (2001), are advancing rapidly and hopefully 
will soon allow us to include explicitly most or all of the physics outlined here. 
Second, some recent insights suggest that the full solution for modified shocks 
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may turn out to be robust after all, and that all of the complications work 
together in a way to find a "critical" solution. 

For example, Malkov (1999) demonstrated recently the existence of a simi- 
larity solution for a steady, strongly modified CR shock. He found a solution to 
the coupled diffusion-convection and Euler's equations in which the flow is highly 
modified, even allowing the total compression to become arbitrarily large. Those 
were exactly the kinds of situations where concerns were raised above about the 
development of non-homologous spectra dependent on many details. However, 
despite the highly nonlinear character of Malkov's similarity solution, the CR 
distribution function takes a simple power-law of the form, f{p) oc inde- 
pendent of the total compression and the form of a given k{p) oc p^" , so long as 
/?s > 0.5. In fact, the flow profile adjusts to the given form of k, so that the par- 
ticle distribution takes this "universal form" . In that work Malkov additionally 
confirmed from the diffusion-convection equation the existence of three distinct 
solutions analogous to those mentioned earlier, giving very different results for 
the efficiency of the acceleration and consequent shock modification. This cal- 
culation also established that no solution with a vanishingly small gas subshock 
can result from the diffusion-convection equation. Malkov identified a critical 
parameter determining the character of the solution and established for 

a given Mach number, M, and CR energy injection efficiency, rj, that the exis- 
tence of test-particle and high efficiency solutions depends on the value of p2. 
Small values lead to the test-particle solution, while large values formally open 
up "intermediate efficiency" and the "high efficiency" solutions. 

As an extension of this insight, Malkov, Diamond and Volk (2000) recently 
argued that modified CR shocks may evolve towards a self- organized critical 
state that finds the combination of critical parameters balancing the energetics 
of the shock at an appropriate equipartition among the components; i.e, an 
"attractor state." The critical solution for a given Mach number depends on 
the self-adjustment of the particle injection rate and the maximum momentum, 
P2- Those in turn are coupled to each other and to the growth, propagation and 
dissipation of the resonant scattering waves and the underlying bulk plasma that 
complete the basic physical system. Those authors suggest that the solution 
would correspond to the shock compression and injection rate that give exactly 
one solution. That, in turn depends on the maximum particle momentum, p2. 
Were p2 to increase above its critical value the total compression would increase, 
weakening the subshock and reducing the injection rate. On the other hand if 
increased escape rates reduced p2 the total compression would be expected to 
decrease, leading to an increase in the injection rate. Both effects would serve 
to return the system to its critical point. 

While these fascinating insights are yet to be confirmed by direct simula- 
tions, they do remind us that DSA is one piece of the full physics of collisionless 
shock formation, and encourages us to look for a unified view of the results. Then 
perhaps we can anticipate a relatively simple outcome to match the empirical re- 
sult that CRs do seem to appear widely with power-law spectra, implying, when 
viewed as a whole, that the details are not so crucial after all. That "whole" 
is currently beyond our understanding, but we are making good strides towards 
clarifying it. 
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